Hylogenesis: 

A Unified Origin for Baryonic Visible Matter and Antibaryonic Dark Matter 
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We present a novel mechanism for generating both the baryon and dark matter densities of the 
Universe. A new Dirac fermion X carrying a conserved baryon number charge couples to the 
Standard Model quarks as well as a GeV-scale hidden sector. CP-violating decays of X, produced 
non-thermally in low-temperature reheating, sequester antibaryon number in the hidden sector, 
thereby leaving a baryon excess in the visible sector. The antibaryonic hidden states are stable dark 
matter. A spectacular signature of this mechanism is the baryon-destroying inelastic scattering of 
dark matter that can annihilate baryons at appreciable rates relevant for nucleon decay searches. 



I. Introduction: Precision cosmological measurements 
indicate that a fraction fif, ~ 0.046 of the energy content 
of the Universe consists of baryonic matter, while £ld — 
0.23 is made up of dark matter (DM) [l|. Unfortunately, 
our present understanding of elementary particles and 
interactions, the Standard Model (SM), cannot account 
for the abundance of either observed component of non- 
rclativistic particles. 

In this Letter we propose a unified mechanism, 
hylogenesis 1 , to generate the baryon asymmetry and the 
dark matter density simultaneously. The SM is extended 
to include a new hidden sector of states with masses near 
a GeV and very weak couplings to the SM. Such sectors 
arise in many well-motivated theories of physics beyond 
the SM, and have received much attention within the 
contexts dark matter models Q, and high luminosity, 
low-energy precision measurements Q . 

The main idea underlying our mechanism is that some 
of the particles in the hidden sector are charged under 
a generalization of the global baryon number (B) sym- 
metry of the SM. This symmetry is not violated by any 
of the relevant interactions in our model. Instead, equal 
and opposite baryon asymmetries are created in the vis- 
ible and hidden sectors, and the Universe has zero total 
B. These asymmetries are generated when (i) the TeV- 
scale states X\ and its antiparticle X\ (carrying equal 
and opposite B charge) are generated non-thermally in 
the early Universe (e.g., during reheating), and (ii) X\ 
decays into visible and hidden baryonic states. The X\ 
decays violate quark baryon number and CP, and oc- 
cur away from equilibrium. Both the visible and hidden 
baryons are stable due to a combination of kinematics 
and symmetries. The relic density of the hidden baryons 
is set by their asymmetry, and they make up the dark 
matter of the Universe. We compute the baryon and 
dark matter densities within a concrete model realizing 



this mechanism in Section II. 

A potentially spectacular signature of our model is that 
rare processes can transfer baryon number from the hid- 
den to the visible sector. Effectively, antibaryonic dark 
matter states can annihilate baryons in the visible sector 
through inelastic scattering. These events mimic nucleon 
decay into a meson and a neutrino, but are distinguish- 
able from standard nucleon decay by the kinematics of 
the meson. In Section III, we discuss this signature in 
more detail, along with its implications for direct detec- 
tion and astrophysical systems. 

We note that our scenario shares some elements with 
Rcfs. (I- 11 1, Du t involves a different production mecha- 
nism and unique phenomcnological consequences. 
II. Genesis of Baryons and DM: In our model, the 
hidden sector consists of two massive Dirac fermions X a 
(a = 1,2, with masses mx 2 > mx 1 > TcV), a Dirac 
fermion Y, and a complex scalar $ (with masses my ~ 
m$ ~ GeV). These fields couple through the "neutron 
portal" (XU C D C D C ) and a Yukawa interaction: 



if D j± X a P R d u c P R d + ( a X a Y c 3>* 



h.c. 



(1) 



Many variations on these operators exist, corresponding 
to different combinations of quark flavors and spinor con- 
tractions. With this set of interactions one can define a 
generalized global baryon number symmetry that is con- 
served, with charges Bx = —(By+B^,) = 1. The proton, 
Y, and $ are stable due to their B and gauge charges if 
their masses satisfy 



I my — m$ I < m p + m e < my + m$ 



(2) 



1 From Greek, hyle "primordial matter" + genesis "origin.' 



Y and $ are the "hidden antibaryons" that comprise the 
dark matter. Furthermore, there exists a physical CP- 
violating phase arg(A* A2C1C2 ) that cannot be removed 
through phase redefinitions of the fields. 

We also introduce a hidden U(l)' gauge symmetry un- 
der which Y and <E> have opposite charges ±e', while X a is 
neutral. We assume this symmetry is spontaneously bro- 
ken at the GeV scale, and has a kinetic mixing with SM 
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FIG. 1: Tree-level and one-loop graphs for decay Xi - 

hypercharge U(l)y via the coupling —^B^Z', 
B 



■> udd. 
where 

^ and Z' are the t/(l)y and t/(l)' field strength ten- 
sors. At energies well below the electroweak scale the 
effect of this mixing is primarily to generate a vector 
coupling of the massive Z' gauge boson to SM particles 
with strength — cw^Qem&- The GeV-scale Z' masses 
we consider here can be consistent with observations for 
10~ 6 <k<1(T 2 0. 

Baryogenesis begins when a non-thermal, CP- 
symmetric population of X\ and X\ is produced in the 
early Universe. These states decay through X\ — > udd 
or Xi — > F$* (and their conjugates). An asymmetry 
between the partial widths for X\ — > udd and X\ — > udd 
arises from interference between the diagrams shown in 
Fig. [I] and is characterized by 

e = — !— [r(Ai udd) - T(X 1 — » udd)] (3) 

256tt 3 | Ci | 2 M i m X2 ' 

where we have assumed that the total decay rate Tx x 
is dominated by X\ — > F$* over the three-quark mode, 
and that mx 2 3> mxi ■ For e ^ 0, Ai decays generate 
a baryon asymmetry in the visible sector, and by CPT 
an equal and opposite baryon asymmetry in the hidden 
sector. These asymmetries can be "frozen in" by the 
weakness of the coupling between both sectors. 

We model the non-thermal production of X\ as a re- 
heating process after a period where the energy content 
of the Universe was dominated by the coherent oscilla- 
tions of a scalar field tp. This field could be the inflaton, 
or it could be a moduli field arising from an underlying 
theory with supersymmetry [l2l ] or a compactification of 
string theory [13J. As tp oscillates, it decays to visible 
and hidden sector states reheating these two sectors. We 
suppose that a fraction of the tp energy density p v is con- 
verted into Ai, X\ states, while the remainder goes into 
visible and hidden sector radiation which quickly ther- 
malizes due to gauge interactions. 

The dynamics of hylogenesis and reheating are gov- 
erned by the Boltzmann equations 

d 



dt 
d_ 

dl 
d_, 

dl 



a 3 p v ) 



v 3 

J- tp a Pip ■ 



(a 3 s) = + T ip a 3 p v /T , 



a 3 n B ) 



iN'xFipO^ptp/m^ 



(4a) 
(4b) 
(4c) 



with tp mass m v and decay rate T v . s = ,shs + ssm = 
(27r 2 /45)g s T 3 is the total entropy density of SM and HS 



states (assumed in kinetic equilibrium at temperature 
T with an effective number of entropy degrees of free- 
dom g s (T)), and ns is the baryon number density in 
the visible sector (i.e. quarks). The scale factor a(t) is 
determined by the Friedmann equation H 2 = (a/a) 2 = 
(8ttG/3) {p v + p r ), where p r = (7r 2 /30).gT 4 is the total 
radiation density and g(T) is the effective number of de- 
grees of freedom. Mx is the average number of X\ states 
produced per tp decay. 

Eq. |4aj) describes the depletion of the oscillating field 
energy due to rcdshifting and direct tp decays and has the 
simple solution p v oc e~ Tvt a~ 3 , while Eq. (|4b|) gives the 
rate of entropy production due to decays and describes 
the reheating of the Universe. We adopt the convention 
that reheating occurs at temperature Trh , defined when 
Pt{Trh) = Plp{Tr,h)- This occurs near the characteristic 
decay time t ~ T" 1 , where the total decay width T v takes 
the form [|, T v = m 3 /(47rA 2 ). Here, A is a large 
energy scale corresponding to the underlying ultraviolet 
dynamics. For example, A - M P1 = 2.43 x 10 18 GeV 
for many moduli in string theory or supergravity. At 
reheating, the radiation temperature is approximately Q 
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100 TeV 



3/2 



(5) 



We require Tru <] 5 MeV to maintain successful nucle- 
osynthesis. 

Eq. (|4cp determines the comoving density of visible 
baryons. The remnant of the intermediate Ai stage ap- 
pears in the right-hand-side of Eq. (|4"c"l). The factor e 
encodes the Ai decay asymmetry. In writing Eq. Q we 
implicitly take mx 1 3> T and T x t S> T v , H. The former 
condition implies inverse decays and scattering reactions 
that could wash out the asymmetry, such as uXi — > dd , 
are suppressed by Boltzmann factors of e~ mXl ^ T ', while 
the latter condition is satisified for |Ci| ^ m ' 2 ,/{ m x 1 ^)- 
The hidden- visible baryon asymmetry can also be washed 
out by F$ — > 3q scattering. A sufficient condition for this 
washout process to be ineffective is 



-1/5 



T RH < (2 GeV) 



^ M 4 mx a m Xb 



(0) 



The allowed Trh increases roughly linearly with the mass 
scale (M 4 m Xl 2 ) 1/6 ■ 

The resulting baryon asymmetry today is given by 



n B /s= j(m v T l p) . (J) 



Assuming that reheating occurs instantaneously, one can 
show analytically that / = 3/4. A numerical solution to 
Eqs. (g]) reveals / ~ 1.2, with less than 10% variation 
over a wide range of (m v ,T v ). Larger values of Trh 
(larger m v for fixed A) allow for greater production of 
baryons. 
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For the parameter values m v = 2000 TcV, A = Mpi, 
N x = 1, we find Trjj ~ 400 MeV and <q B /e ^ 2.5 x 10~ 7 . 
The observed value of the baryon asymmetry is obtained 
for lm[\l\ 2 C 2 /(i}m Xi / (M 4 m X2 ) ~ 3. Smaller values of 
e and m v are viable for A < Mj>\. 

We have implicitly assumed that the Z' maintains ki- 
netic equilibrium between the SM and hidden sectors. 
This will occur if Tz'(l/'y) > H, where 7 is a relativistic 



time dilation factor, which implies [14 1 



« > 1.5 x 10" 8 



1/2 



m z > 
GeV 



T 
GeV 



3/2 



(8) 



provided Trh > After baryogenesis, the CP- 

symmetric densities of hidden states are depleted very 
efficiently through annihilation YY — > Z'Z 1 and $$* — > 
Z'Z' provided mz' < my, m$, with the Z' decaying 
later to SM states by mixing with the photon. The cross- 
section for YY ->• Z'Z 1 is given by [3] 



(crv} 



16tt m Y 



1 — m%,/m Y 



(9) 



(1.6 x 10- 25 cm d /s) 



0.05 



/ 3GeV\ 2 
my J 



Annihilation of $*$ is given by a similar expression. 
These cross sections are much larger than what is needed 
to obtain the correct DM abundance by ordinary ther- 
mal freeze-out, and all of the non-asymmetric DM den- 
sity will be eliminated up to an exponentially small re- 
mainder (l5j |. Note that the annihilation process may 
occur later than is typical for thermal freeze-out for 
Trh ^ TOy,$/20, but even in this case the remaining 
non-asymmetric density will be negligibly small (l6| . 

The role of the hidden Z' in our model is to ensure 
the thermalization and symmetric annihilation of Y and 
$. A more minimal alternative is to couple $ to the SM 
Higgs boson h via the operator £ |/i| 2 |<I>| 2 . For £ > 10~ 3 , 
this interaction, together with FX$ and |<&| 4 , appears 
to be sufficient for both thermalization and symmetric 
annihilation. 

The residual CP-asymmetric density of Y, $ is not 
eliminated and makes up the DM [ItJ. The relic num- 
ber density is fixed by total baryon number conservation: 
iiy = n$ = Tig. Thus the ratio between the energy den- 
sities of DM and visible baryons is 



Qd/Qb = {my + m$)/m p 



(10) 



Present cosmological observations imply fid /fib = 
4.97 ± 0.28 [lj, which corresponds to a range 
4.4 GeV < my +771$ < 4.9 GeV, or 1.7 GeV < 
my,m$ < 2.9 GeV when combined with the con- 
straint I my — m$| < m p + m e . 

III. Dark Matter Signatures: A novel signature 
of this mechanism is that DM can annihilate nuclc- 
ons through inelastic scattering processes of the form 



P 



K+ 



FIG. 2: Diagram for induced nucleon decay processes pY — ► 
and ->■ K + Y. 

YN ->• $*M and $V -> YM mediated by X lj2 , where N 
is a nucleon and M is a meson (Fig. [2]) . We call this pro- 
cess induced nucleon decay (IND). IND mimics standard 
nucleon decay N — > Mv, but with different kinemat- 
ics of the daughter meson, summarized in Table HI For 
down-scattering processes, where the mass of the initial 
DM state is greater than the final DM state, the meson 
momentum pm from IND can be much greater than in 
nucleon decay. The quoted range of pm corresponds to 
the range of allowed masses (my,m$) consistent with 
Eqs. (IS [10]). For fixed masses, pm is monochromatic, 
with negligible broadening from the local DM halo ve- 
locity. (We also note a related study considering lepton- 
number violating inelastic DM- nucleon scattering (18|.) 

To estimate the rate of IND we consider the specific op- 
erator (X a /M 2 )(u c d)ji(Xs)R that mediates <&p -)• YK + , 
illustrated in Fig. [2] Treating the $ and Y states as 
spectators, the hadronic matrix element can be estimated 
from the value computed for the p — > K + v decay through 
the corresponding three-quark operator fioj ] . We find 
that the sum of the IND scattering rates p<I> — > K + Y 
and pY K + Q* is given by 



(crv) IND = C(10" 39 cm7s) 



TeV d 



try 



m a M 2 /\* a Ca 



(11) 



where 0.5 < C < 1.6, depending on m$.y within the al- 
lowed range. We expect IND modes from other operators 
to be roughly comparable. This estimate, which relies on 
a chiral perturbation theory expansion that is expected 
to be poorly convergent for pu ~ 1 GeV, is approximate 
at best. 

An effective proton lifetime t p can be defined as the 
inverse IND scattering rate per target nucleon, t~ x = 
nDAi(o"v)iND- With a local DM density of 0.3 GeV/cm 3 , 
(pv)ind = 10 _39 cm 3 /s corresponds to a lifetime of 
T, p ~ 10 32 yr. This is similar to the current lifetime bound 



Decay mode 


pir (MeV) 


Pm 


(MeV) 


TV ->• 7T 


460 


800 


- 1400 


N -> K 


340 


680 


- 1360 


N -> 7? 


310 


650 


- 1340 



TABLE I: Daughter meson M £ {it, K, 77} momentum pm for 
standard nucleon decay (SND) and down-scattering IND. 
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on p — > K + v of 2.3 x 10 33 yr [20]. However, existing nu- 
cleon decay bounds may not directly apply to IND due to 
the non-standard meson kinematics [21], and additional 
suppression can arise from the second factor in Eq. ([lip . 

There is also a direct detection signal in our model due 
to the hidden Z': Y and <I> can scatter elastically off pro- 
tons. The effective scattering cross-section per nucleon 
for either Y or $ is spin-independent and given by 
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a^ 1 = (5 x 10 



- 39 cm 2 ) 



2Z 
~A 



ao5y Go- 5 ) ( 



GeV/ 
0.1 GeV 

mz< 



(12) 



where ^at is the DM-nucleon reduced mass. For a DM 
mass of 2.9 GeV, this is slightly below the current best 
limit from CRESST j22[. The effective nucleon cross- 
section will be much lower if the hidden vector is replaced 
by a <i>-Higgs coupling. 

Annihilation of DM can generate energetic particles 
that destroy the products of nucleosynthesis [23[, create 



a neutrino flux |24[ , and modify the properties of the cos- 
mic microwave background [25|. In our scenario there is 
almost no direct DM annihilation after freeze-out, but 
similar effects can arise from IND scattering. However, 
the rate for this scattering given in Eq. (fTTj) is many 
orders of magnitude less than the cor resp onding limits 
on DM annihilation cross-sections [23, 25 1, and thus we 
expect that IND will have no visible effect in the cosmo- 
logical setting. 

The effects of IND can become important in astrophys- 
ical systems with very large densities of nucleons such as 
neutron stars and white dwarfs. In both cases, we find 
that the rate of IND in the stellar core typically becomes 
large enough that it reaches a steady state with the rate 
of DM capture through elastic scattering [2(| [2?} ■ This 
has the effect of heating the stellar interior in much the 
same way as DM annihilation. Nucleons are also de- 
stroyed, but the number is only a tiny fraction of the to- 
tal. Current observations of white dwarfs constrain stel- 
lar heating by DM or IND, but the bounds depend on the 
local DM densities which are not known precisely 2^, 2$ ■ 
These bounds are much weaker than for monopole catal- 
ysis of nucleon decay since the anti-DM product of 
IND is unable to destroy any more nucleons. 
IV. Conclusions: We have presented a novel mecha- 
nism to generate dark matter and baryon densities si- 
multaneously. Decays of a massive X\ state split baryon 
number between SM quarks and antibaryons in a hid- 
den sector. These hidden antibaryons constitute the dark 
matter. An important signature of this mechanism is the 
destruction of baryons by the scattering of hidden dark 
matter. 
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